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1.6 Relations between stress and rate-of-strain tensors 

When the fluid is at rest on a macroscopic scale, no tangential stress acts on a surface. There 
is only the normal stress, i.e., the pressure — pSij which is thermodynamic in origin, and is 
maintained by molecular collisions. Denoting the additional stress by which is due to the 
relative motion on the continuum scale, 

<7ij = -pSij + Tij (1.6.1) 

The second part is called the viscous stress and must depend on gradients of velocity, 
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1.6.1 Newtonian fluid 



For many fluids in nature such as air and water, the relation between and dqi/dxi are 
linear under most circumstances. Such fluids are called Newtonian. For second-rank tensors, 
the most general linear relation is, 

Tij = Cijl m ^-. (1.6.2) 

<JX m 

where Cije m is a coefficient tensor of rank 4. In principle there are 3 4 = 81 coefficients. It 
can be shown (Spain, Cartesian Tensors) that in an isotropic fluid the fourth-rank tensor is 
of the following form: 

Cijim = ^ijSi m + ^(SuSjm + 5 im 5j£) (1.6.3) 

Eighty one coefficients in CV^ m reduce to two: A and /i, and 

T -=^{dx- + dx-J +X dx- e ^ (1 - 6 - 4) 

where /i, A are viscosity coefficients depending empirically on temperature. 
Note that the velocity gradient is made up of two parts 

dqj = 1 / dqj | dqA | 1 / % dqA 
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where 



is the rate of strain tensor, and 



is the vorticity tensor. Note also that (1.6.4) depends only on the rate of strain but not on 
vorticity. This is reasonable since a fluid in rigid-body rotation should not experience any 
viscous stress. In a rigid-body rotation with angular velocity u>, the fluid velocity is 



q = to x r 
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The vorticity components are not zero; for example, 

2^12 = TT^- ~ TT^ = TT- (^2^3 ~ ~ (^3^1 - W1X3) = ~ 2uJ S 

OX 2 OXi OX 2 OXi 

Hence r^- cannot depends on Vt^ and only on e^. 



The trace of is 



a ti = (2/i + 3A)J^ = (3A + 2/x)V • q. 



where k = 3A + 2/x = is called the bulk viscosity. 

For incompressible fluids V • q = 0; the viscous stress tensor is 



( dqi dqA 

a «=»[dx- + dx- i ) (L6 - 7) 

The total stress tensor is therefore 

The governing equations for an incompressible Newtonian fluid may now be summarized: 

— p = 0, (incompressibility) (1.6.9) 

|^ = 0, (continuity) (1.6.10) 

Dqi dp d 2 q,i 

p—— = — 7: h a— — h pfi, (momentum conservation) (1.6.11) 

Dt OXi oxjOXj 
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after using continuity. The last equation (1.6.11) and sometimes the set of equations (1.6.11), 
(1.6.9) and (1.6.10), is called the Navier-Stokes equation(s). Now we have just five scalar 
equations for five unknowns p,p, and g«. Boundary and initial conditions must be further 
specified. For example on the surface of a stationary rigid body, no slippage is allowed, so 
that 

qi = 0, on a rigid stationary surface (1.6.12) 



1.6.2 Non-Newtonian fluids 



Many fluids such as toothpaste, gel, honey, heavy oil (DNAPL), etc., flows like a fluid if the 
shear stress is above a critical value, and behaves like a solid if below. Of geological interest 
is the mud which is a mixture of water with highly cohesive clay particles. From volcanic 
eruption, lava can mix with rain, melting snow, or lake water to form mud, which flows down 
the hill slope, carries along stones, trees and other debris, to cause severe damages. In some 
mountainous areas, heavy rainfall infiltrates the top soil and causes mud to slide and flow 
into rivers to form hyper concentrated fluid-mud. These fluids are called non-Newtonian 
since the relation between and is nonlinear. 

We shall only discuss a special model of non- Newtonian fluid, i.e., the Bingham plastic 
model. For simple shearing flow u = u(y), the constitutive relation for a Bingham plastic is 

du 
dy 

du I 



= 0, r < r c ; 



(1.6.13) 
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T > T, 



where r c is called the yield stress and \x the Bingham viscosity, both of which depend on the 
clay concentration C . 

In three dimensions, the Bingham model can be generalized by introducing the second 
invariants of the stress and rate-of-strain tensors. The second invariant of the viscous stress 
tensor is 
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(1.6.14) 



2.3 



'31 



(T11T22 + T11T33 + T22T33) 



Similarly the second invariant of the rate of strain tensor is 
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(1.6.15) 
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The Bingham plastic law is then 

0, if ^flh<T c , (1.6.16) 



T i:j = 2fieij + T c -?=i!=, if \[u~ T > t c . 
V He 

This is due to Hohenemser and Prager (1936). 

In simple shear u = u(y), the only non-zero components of r^- and are r xy and e xy . 
The Bingham law reduces to 

e xy = 0, if \r xy \ < t c , (1.6.17) 



In other words, 

du 



. = 0, if \r xy \ <r c , (1.6.18) 

du ( du\ . 

= A«^- + r c sgn — , if |rxj,| > r c . 



ay \oy 
We shall examine some examples later 



